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ABSTRACT 

Wavelet analysis can be used to measure the power spectrum of solar wind fluctuations along a line in any 
direction (6, 4>) with respect to the local mean magnetic field Bq. This technique is applied to study solar wind 
turbulence in high-speed streams in the ecliptic plane near solar minimum using magnetic field measurements 
with a cadence of eight vectors per second. The analysis of nine high-speed streams shows that the reduced 
spectrum of magnetic field fluctuations (trace power) is approximately azimuthally symmetric about B in 
both the inertial range and dissipation range; in the inertial range the spectra are characterized by a power-law 
exponent that changes continuously from 1.6 ± 0.1 in the direction perpendicular to the mean field to 2.0 ± 0.1 
in the direction parallel to the mean field. The large uncertainties suggest that the perpendicular power-law 
indices 3/2 and 5/3 are both consistent with the data. The results are similar to those found by Horbury et 
al. (2008) at high heliographic latitudes. Comparisons between solar wind observations and the theories of 
strong incompressible MHD turbulence developed by Goldreich & Sridhar (1995) and Boldyrev (2006) are not 
rigorously justified because these theories only apply to turbulence with vanishing cross-helicity although the 
normalized cross-helicity of solar wind turbulence is not negligible. Assuming these theories can be generalized 
in such a way that the 3D wavevector spectra have similar functional forms when the cross-helicity is nonzero, 
then for the interval of Ulysses data analyzed by Horbury et al. (2008) the ratio of the spectra perpendicular and 
parallel to Bq is more consistent with the Goldreich & Sridhar scaling P±/Pu oc v 1 ^ than with the Boldyrev 

scaling v 1 ! 2 . The analysis of high speed streams in the ecliptic plane does not yield a reliable measurement of 
this scaling law. The transition from a turbulent MHD-scale energy cascade to a kinetic Alfven wave (KAW) 
cascade occurs when k±pi ~ 1 which coincides with the spectral break. At slightly higher wavenumbers, in 
the dissipation range, there is a peak in the power ratio with P±/P\\ ^> 1. The decay of this peak may be 
caused by the damping of KAWs which is predicted to occur near k±pi ~ 4. 
Subject headings: Solar wind — turbulence, magnetohydrodynamics, scaling laws 



1. INTRODUCTION 

In the presence of a strong mean magnetic field Bq, MHD 
turbulence is spatially anisotropic: velocity and magnetic field 
fluctuations vary more rapidly in the direction perpendicular 
to Bo than parallel to Bo. It is important to recognize that 
fluctuations at a given scale are most strongly affected by the 
local mean magnetic field at roughly 3 to 5 times that scale. 
When this scale dependence of the local mean magnetic field 
is taken into account, the anisotropy of the turbulence be- 
comes scale dependent because the direction of the local mean 
magnetic field changes with scale. Thus, anisotropy is a local 
property of the turbulence that depends on both position and 
scale. 

In light of this, it is natural to expect wavelet analysis to be 
an effective tool for the study of anisotropy in plasma turbu- 
lence. Wavelet analysis allows a signal fit") to be decomposed 
into components that are localized in both time and frequency 
(or wavelet scale). Therefore, unlike the Fourier transform, 
the wavelet t ransform is well adapted to the study of tran- 
sient signals (Daubechies 1992). Because of these properties, 
wavelet analysis can be used to study the power spectrum of 
turbulent fluctuations as a function of th e direction of the loca l 
mean magnetic field as first shown by Horburv e t al.l d2008l) 
using solar wind data. 

Horbury et al. (2008) studied a 30 day interval of 1 sec- 
ond magnetic field data acquired by the Ulysses spacecraft at 
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high heliographic latitudes during solar minimum conditions 
in 1995 at a heliocentric distance of ~ 1.4 AU. The high lati- 
tude flow of the solar wind at this time can be characterized as 
high speed, highly Alfvenic (large magnitude of the normal- 
ized cross-helicity), with the average magnetic field directed 
roughly radially outward at the largest scales. By binning the 
square amplitude of the wavelet coefficients at a given scale 
according to the direction of the local mean magnetic field 
Bq at that scale they estimated the average power as a func- 
tion of frequency and of the direction of Br, . They found that 
the power law index of inertial range fluctuations varies from 
approximately 5/3 to approximately 2 as the angle 9 between 
Bq and the average flow direction decreases from ir/2 to 0. 

This was interpreted by Horbury et al. (2008) as evidence 
that the 3D wavevector s pectrum of solar wind fluct uations 
has the form given in the iGoldreich & Sridhar ( 1995) theory 
of strong incompressible MHD turbulence. This interpreta- 
tion is somewhat paradoxical, however, because the Goldre- 
ich & Sridhar (1995) theory does not apply to turbulence with 
non-vanishing cross-helicity as is usually found in the solar 
wind. Nevertheless, the results are intriguing because they 
demonstrate that the power law index of solar wind fluctua- 
tions varies with the angle 9 between the flow direction and 
the local mean magnetic field, at least for the one data record 
studied by Horbury et al. (2008). This has never been seen 
before. In fact, based on earlier work bv lSari & Vallevl ( 1 19761) 
and possibly others, it is a common belief that the spectral ex- 
ponent of solar wind fluctuations is independent of the angle 
9; a view confirmed in a recent study by Tessein et al. (2009). 
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However, the study by Tessein et al. (2009) does not take 
into account the scale dependence of the local mean magnetic 
field and this may explain, in part, why no dependence of the 
spectral exponents on was found in that study. 

Near solar minimum, the solar wind in the ecliptic plane is 
composed of recurring high-speed streams interspersed with 
lower speed wind. These streams originate in equatorial coro- 
nal hole s which extend from the pola r regions of the sun to the 
equator dZirkerll 977; Schwenn 2006) and have physical char- 
acteristics that are similar to the high-speed, high-latitude, so- 
lar wind found above polar coronal holes near solar minimum. 
Therefore, the variation in the power law index seen in high 
latitude Ulysses data by Horbury et al. (2008) should also 
be present in high-speed streams in the ecliptic plane around 
solar minimum. The purpose of this study is to investigate 
this expectation using magnetometer data from NASA's two 
Stereo spacecraft and to interpret the results using anisotropic 
theories of incompressible MHD turbulence. It should be em- 
phasized that the physical quantity measured in this study is 
the so called reduced spectrum which is an integral of the 3D 
wavevector spec trum over a plane perpendicular to the mea n 
magnetic field dMatthaeus & Goldstein|[T982l; lPodestall2009l) . 

2. DATA SELECTION 

NASA's two identically equiped Stereo spacecraft, Stereo 
Ahead (STA) and Stereo Behind (STB), were launched on 
25 October 2006. The spacecraft both follow an earth- 
like orbit with one spacecraft moving progressively ahead 
of the earth and the other moving progressively behind the 
earth. This study uses dat a from the Stereo flux gate magne- 
tometer (lAcufia et alJl2008h part of the IMPACT instrument 
dLuhmann et alj|2008l) . The data has a uniform cadence of 8 
vectors per second and contains no data gaps or missing data 
over the time intervals studied here. 

Between April 2007 and June 2008 the solar wind in the 
ecliptic plane consisted a succession of high-speed streams 
of varying amplitudes interspersed with low-speed wind — 
typical conditions near solar minimum. The data used in 
this study is restricted to high-speed streams; the leading 
and trailing edges of high-speed streams were usually ex- 
cluded. High-speed streams of long duration with sustained 
high speeds were prefered over shorter lived more rapidly de- 
caying streams because they provided larger contiguous sam- 
ples of high speed wind. Streams were identified using plasma 
data, proton velocity and density, from t he PLASTIC instru - 
ment on board the two Stereo spacecraft dGalvin et al.l l2008). 
The high-speed streams analyzed in this study are displayed in 
Figure [T] and their properties are listed in Table 1. The length 
of each data record is typically 3 or 4 days. 

Near solar minimum, high speed streams are produced by 
equatorial coronal holes and are believed to be similar to 
the high speed wind at high heliographic latitudes. A sin- 
gle stream is usually embedded in a single magnetic sector of 
the interplanetary magnetic field and contains no sector cross- 
ings. For the nine intervals listed in Table 1, this was checked 
using running averages of the magnetic field vector to cre- 
ate 2D scatter plots of the average values of (Br, Bt) versus 
time. Averaging windows with durations of 20, 40, and 60 
minutes were used and, for simplicity, the offset from the be- 
ginning of one averaging window to the beginning of the next 
window was half the window width. The data is assumed to 
be contained in a single magnetic sector if the scatter plot re- 
mains in the same quadrant, either the second (Bq inward) or 



the fourth (Bq outward), for a majority of the time with no 
points or very few points in the opposite quadrant (inward or 
outward). 

The only other criteria used in the data selection is the re- 
quirement that power spectra obtained by the wavelet anal- 
ysis described below should form well defined power laws 
throughout the inertial range and at all angles of interest. 
In one instance, some of the spectra obtained using data 
from December 2007 showed unusual fluctuations and erratic 
points suggesting there may have been inadequate sampling 
of the power in some angle bins or unusual statistical fluctu- 
ations in the time series, possibly caused by solar wind tran- 
sients. Therefore, this interval was discarded. 

3. WAVELET ANALYSIS TECHNIQUE 

3.1. Continuous wavelet transform 

The continuous wavelet transform of a function f(t) is 
defined by 

' T -t' 



F(s,t) 



-1/2 



f(r)dr, (1) 



where s is the wavelet scale, t is the time, ip(t) is the 
mother wavelet, /_ KKOI 2 dt = 1, and the asterisk '*' de- 
notes the complex conjugate. The wavelet scale s is related 
to the Fourier frequency as discussed below. The continu- 
ous wavelet transform allows any square integrable function 
f(t) to be written as a continuous superposition of scaled 
and time-shifted wavelets by means of the inversion formula 
(Daubechies 1992) 
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(4) 



is the Fourier transform of tp(t). The distribution of energy in 
the (s, t) plane (related to the time-frequency plane) is indi- 
cated by the relation 

r w)?dt=±r rV(-,t)i 9 ^, (5) 

J — oo J— oo J — oo 

which is similar to Parseval's formula. 

In the special case where f(t) is real valued and ip(cu) van- 
ishes for ui < 0, the inversion formula takes the simplified 
form 



,^ 2 f°° ds 



dr Re 



-1/2 



F(s,t) 



(6) 

where C is given by ((3). Note that only positive scales s > 
are included in the integral ©. Parseval's relation © takes 
the form 

/OO c\ rOO 7 /"OO 

\f(t)\ 2 dt=^ J/ dt\F(s,t)\ 2 . (7) 
-oo u JO S J -oo 

Equations (0 and © are derived in Appendix A. 
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FIG. 1 . — Speed profiles obtained from the plasma instrument (PLASTIC) on Stereo. The time intervals used in the wavelet analysis lie between the two vertical 
lines. The labels 1 through 9 correspond to the time intervals in Table 1. 



TABLE 1 

High-speed streams in the ecliptic plane near solar minimum 



# 


Year 


Begin 


End 


Days 


SC a 


V! (Hz) b 




VI (Hz) b 




Perpendicular 
Exponent 


Exponent 
-All data d 


1 


2007 


28 Apr 00:00 


01 May 00:00 


3 


STB 


1 x 10" 


-_> 


2 x 10" 


i 


1.60 ±0.07 


1.62 ±0.015 


2 


2007 


25 May 00:00 


28 May 02:39 


3.11 


STB 


4 x 10" 


-3 


1 x 10" 


l 


1.55 ±0.07 


1.58 ±0.015 


3 


2007 


27 Aug 12:00 


30 Aug 12:00 


3 


STB 


7 x 10" 


-3 


3 x 10" 


l 


1.70 ±0.07 


1.63 ±0.015 


4 


2007 


15 Nov 00:00 


18 Nov 00:00 


3 


STA 


7 x 10" 


-3 


2 x 10" 


l 


1.55 ±0.07 


1.61 ±0.015 


5 


2008 


08 Jan 00:00 


11 Jan 00:00 


3 


STA 


5 x 10" 


-3 


2 x 10" 


l 


1.52 ±0.07 


1.57 ±0.015 


6 


2008 


13 Feb 00:00 


18 Feb 00:00 


5 


STA 


7 x 10" 


-3 


2 x 10" 


l 


1.65 ±0.07 


1.60 ±0.015 


7 


2008 


08 Mai' 00:00 


1 1 Mai' 00:00 


3 


STB 


1 x 10" 


-2 


2 x 10" 


l 


1.55 ±0.07 


1.61 ±0.015 


8 


2008 


04 Apr 00:00 


08 Apr 00:00 


4 


STB 


5 x 10" 


■3 


2 x 10" 


l 


1.51 ±0.07 


1.62 ±0.015 


9 


2008 


02 May 12:00 


06 May 00:00 


3.5 


STB 


7 x 10" 


■3 


2.5 x 10 


-i 


1.61 ±0.07 


1.63 ±0.015 



a Spacecraft: STA = Stereo A, STB = Stereo B 

b The fits used to determine the power law exponents as a function of the angle 9 were performed over the frequency interval 

v\ < v < V2- 

c Power law index when 9 ~ 90 degrees obtained from the results in Fig. [6] 

d Power law index over the frequency interval 10~ 3 < v < 10 _1 Hz obtained using all data for the indicated time interval. 
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The mother wavelet used in this study is the Morlet wavelet and the inversion formula © can be approximated by 



1p(t) = 7T 



_ „ — l/4r itv t 



» 2 /2 ]p -t 2 /2 



where, in this study, u>q = 6. Its Fourier transform is 



lj,(cu) = 2 1 /2^1/4 [e -(— o) 2 /2 _ e -^/2 f 



2 /2l 



(8) 



(9) 



M-l JV-1 



To leading order, the Morlet wavelet (0 is a Gaussian en- 
velope modulated by a complex exponential. This form is 
chosen because it is well localized in both time and frequency 
with a small time-frequency uncertainty product. If small cor- 
rection terms are introduced to make i)(u>) vanish for u> < 0, 
then the inversion formula © may be used. The magnitude 
of these correction terms are small and can often be neglected 
in practice. For the Morlet wavelet ((8), numerical evaluation 
of the constant C in equation (01 yields the value C 1.06 
when ojo = 6. 

To relate the wavelet scale s of the Morlet wavelet to an 
equivalent Fourier frequency u>, recall that the mean fre- 
quency u> of the function ij}{uj) (assumed positive valued) sat- 
isfies 



(uj — u)tp(w) du> = 0. 



(10) 



Substituting the function (O into this equation shows that 
w ~ wo is an accurate approximation when u>o ^> 1. Because 
the Fourier transform of the scaled wavelet is proportional to 
ip(soj), it follows that the wavelet scale s corresponds to the 
mean Fourier frequency 



2ns 7 



(11) 



where v is the frequency in Hz and lo = 2ttv is the frequency 
in radians per second. A more accurate approximation is ob- 
tained by computing the continuous wavelet transform of the 
function f(t) = cos(ojt) and then finding the scale s where 
\F(s, t)\ 2 is a maximum. To leading order, this yields the re- 
fined approximation derived in appendix B 



2ns 



1 



1 



(12) 



3.2. Practical implementation 

Suppose the function f(t) vanishes outside the interval 
< t < T. If f(t) is sampled at the discrete times t n = nAt, 
where At = T/N and n = 0,1,2, ... ,N - 1, then the 
wavelet transform (Q]) can be approximated by the Riemann 



F(s,kAt) 



N-l 

E s- 1 ? 

71=0 



(n - k)At 



f {nAt) At. (13) 



This assumes that At is small enough to resolve the most 
rapid variations of the function f(t). Likewise, the wavelet 
functions must be sufficiently well resolved at all scales s 
of interest (at least four samples per wave period). Let 
s m = Soa m , where so is the smallest scale, a > 1, and 
m = 0,1,2, ... ,M — 1. Then As = s m = (a~l)s m 



f(nAt) - ^ E E Re i s ^ 



m=0 k=0 



(n - k)At 



x F(s m ,kAt) 



] («-!) 



At 



(14) 



For this approximation to be valid the minimum and max- 
imum values of s must include the smallest and largest 
timescales of the signal f(t n ) and the grid size M along the 
s axis must be chosen large enough to provide adequate fre- 
quency resolution. To resolve the scaled wavelet function in 
equation ( fT3l ) requires at least four samples per wavelet period 
at the minimum scale sq. This implies z/ max < 1/(4 At) and, 
therefore, using equation ( fTTT ), sq > 2ujoAt/TT ~ AAt. The 
maximum scale is typically less than the record length T. 

From equation (0, the total energy of the signal is ex- 
pressed in terms of the wavelet coefficients by the Parseval- 
like relation 



N-l 

E 

n=0 



M-l N-l 



\f( tn )fAt^-j2 Ei^w«)i 



,(a- l)At 



m— n— 



(15) 

The right-hand side of equation ( TT3T > describes the distribution 
of signal energy in the time-frequency plane. In this study, the 
minimum scale is sq — 2uooAt/Ti ~ AAt and the maximum 
scale s max is much less than the record length T in which case 
the right-hand side of ( TTBT l only includes contributions to the 
energy from scales less than or equal to s max - 



The component of th e sig nal at scale s r 
tion to the total energy (TT51) given by 



makes a contribu- 



2 Ar_1 

- Y, \F{s m ,t n ) 
n=0 



, (a - l)At 



(16) 



P{Vm) 



The power spectral density P(v) is defined such that P(v)dv 
is the power in the frequency interval from v to v + dv. Using 
the relation dv ~ ui ds/2irs 2 and the fact that the average 
pow er is equal to the energy divided by the time, equation 
i fToT l implies 

:^£|F( Sm ,t„)| 2 , (17) 
71=0 

where the frequency v m is defined by equation ( TTZb and T is 
the record length. This expresses the Fourier power spectrum 
(frequency spectrum) in terms of the wavelet coefficients. For 
the vector field B(t), the total power is the sum of the power 
of the three orthogonal components so that, in equation ( fTTl ), 
|F(M)| 2 = \F R (s,t)\ 2 + |F T ( S ,t)| 2 + \F N (s,t)\ 2 . 

In practice, the wavelet coefficients F(s m , t n ) may be com- 
puted using the FFT The right-hand side of equation ( fT3l con- 
sists of a convolution of the sequence /„ = f(t n ) with the 
sequence 

g n = i>*{-t n /s). (18) 

Note the minus sign in the argument of the function ^*{t). 
Recall that for periodic sequences of length L, the discrete 
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Fourier transform of the convolution 

L-l 

h n = fn * 9n = ^ fm9n—mi 71. = 0, 1, 2, . . . , L — 1, 

m=0 

(19) 

is equal to iJ^ = FkGk (Briggs & Henson 1995), where 

L-l 

F k = fn exp(~i2nkn/L), k = 0, 1, 2, . . . , L - 1, 

(20) 

is the discrete Fourier transform of the sequence /„ and the 
inverse transform is defined by 

1 

/ n =-V F fc exp(+i2vrfcn/ J L), n = 0, 1, 2, . . . , L-l. 

k=0 

(21) 

The integer L > N must be chosen to avoid aliasing (wrap 
around effects) at the largest wavelet scale s and is typically 
a power of 2 to speed up the FFT algorithm. The envelope of 
the scaled Morlet wavelet exp(— t 2 /2s 2 ) has an approximate 
temporal duration of 6s. Therefore, the required duration for 
the zero padding is 3s max and L must be chosen such that 
(L — N)At > 3s max . Note that the sequence /„ is padded 
with zeros to length L, but g n is not since g n is defined by 
equation ( TT~8T > for — L/2 < n < L/2 and for any length L. 
The wavelet coefficients ( [T3l > at a given scale s are computed 
as follows. Compute the FFT of /„ and g n at that scale, form 
the product FkGk, inverse transform the sequence FkGk, and 
then multiply by s~ x / 2 At. The first N terms of the resulting 
sequence of length L yield the coefficients F(s, t n ) for n = 
0,1,2,. ..,N-l. 

The implementation of the wavelet transform has been 
tested using synthetic data. The wavelet transform of a 1 
Hz cosine wave f(t) = cos(wi) is shown in Figure [2] The 
upper panel shows the magnitude of the wavelet transform 
\F(s m ,t n )\ 2 obtained from the sequence /„ = f(nAt) us- 
ing the sampling time At = 1/200 sec and a record length of 
20 sec. The lower panel shows the exact solution \F(s,t)\ 2 
obtained by evaluating the integral (HJ in closed form (Ap- 
pendix B). Except for the transient near t = 0, the results for 
the wavelet transform agree with the exact solution. For this 
particular calculation the minimum scale is so ~ 8At, the 
maximum scale is s = 5 sec, and M = 200. The equiva- 
lent frequency range extends from v = 0.19 Hz to v = 25 
Hz. The cosine wave is a useful test case and the frequency of 
the wave is easily changed to investigate the behavior of the 
numerical solutions as a function of frequency. 

Many other tests were performed. One more example shall 
now be described. Consider a 20 sec signal f(t) composed 
of a sum of 91 cosine waves having unit amplitude and ran- 
dom phases. The frequencies of the waves extend from 1 Hz 
to 10 Hz in increments of 0.1 Hz. The power spectrum of 
fit) is computed using both the wavelet transform ([T7| > and 
the FFT which is used to compute a frequency spectrum or 
periodogram (smoothed or unsmoothed). The results in Fig- 
ure |3] demonstrate the reasonable agreement between these 
two independent techniques. For the example shown in Fig- 
ure |3] , the wavelet transform was computed using the mini- 
mum scale so ~ 4At, the maximum scale s = 5 sec, and 
M = 200. The FFT spectrum shown in Figure [3] was ob- 
tained using a Papoulis smoothing window with a bandwidth 




2 4 6 8 10 
Time (s) 




2 4 6 8 10 
Time (s) 



FIG. 2. — The wavelet transform of f(t) = cos(wt) obtained using a Mor- 
let wavelet with u>q = 6. The numerical results in the upper plot agree with 
the exact solution in the lower plot. The transient near t = (upper plot) is 
an artifact of the numerical method. 




Frequency (Hz) 



FIG. 3. — The power spectra density (PSD) obtained using the wavelet 
transform (blue) and the FFT (red) for a signal f(t) composed of 91 cosine 
waves having unit amplitude, random phases, and frequencies from 1 Hz to 
10 Hz in increments of 0.1 Hz. The PSD equals 5 Hz -1 between 1 Hz and 
10 Hz (dashed line). 

of 0.1 Hz dPercival & Walden|[l993h . 

4. ANALYSIS OF SOLAR WIND DATA 

The primary goal of this study is to investigate how power 
is distributed as a function of both frequency v and the direc- 
tion (6, cf>) of the local mean magnetic field. The local mean 
magnetic field at time t n and wavelet scale s is obtained by 
weighting the time series with a Gaussian curve centered at 
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time t n . Thus, the mean magnetic field at time t n and wavelet 
scale s is proportional to 



N-l 



B n (s) = ^2 B m exp 



m=0 



2A 2 s 2 



(22) 



where B n = B(t n ). The width of the scaled Gaussian 
curve as measured by its full width at 1/e of its maximum 
is 2 3 / 2 As ~ 2.8As. The dimensionless parameter A deter- 
mines the timescale of the average and in this study A = 1. 
The local mean magnetic field d22l defines a direction (/x, <j>) n 
at each time step n, where /i = cos(0). The average velocity 
of the solar wind lies approximately in the heliocentric radial 
direction defined as 9 — or, equivalently, fi = 1. 

Note that the local mean magnetic field in the work of Hor- 
bury et al. (2008) does not contain the factor of 2 which 
appears in the exponent in equation d22l . The duration or 
timescale of the average (l22l is therefore slightly larger than 
that employed by Horbury et al. (2008). This extraneous fac- 
tor is not expected to significantly change the results of the 
analysis. 

The angular distribution of power may be defined in differ- 
ent ways depending on the physical quantity of interest. The 
power spectral density may be defined as a function of fre- 
quency and angle such that the integral over all solid angles 
and all frequencies is equal to the total power of the signal 
(energy/time). However, in the solar wind, the direction of 
the local mean magnetic field is unevenly distributed over the 
surface of the unit sphere. Consequently, for a given direction 
Bo, the energy at a given frequency is the product of the aver- 
age energy of the fluctuations at that frequency when the field 
points in a solid angle dtt about the direction Bq multiplied 
by the number of times the field Bq points in that direction. 
To describe the distribution of energy of the fluctuations with- 
out the weighting caused by the uneven distribution of direc- 
tions of the local mean magnetic field, it is necessary to use 
the average energy of the fluctuations at that frequency when 
the field points in a given direction (a given angle bin). This 
is the approach adopted by Horbury et al. (2008) and also the 
approach adopted in this study. These two distinct approaches 
are now described in detail. 

4.1. Angular distribution of power 

The power spectral density (T% may be decomposed into 
an angular distribution P(v, [X, 4>) such that the integral over 
the unit sphere yields the power spectral density at frequency 
v. In discrete form this is written 

JV^-lAV-l N-l 

E E p( v m,nk,<i>i) A ^fe = Cuj T E i F ( s ™>v 

(23) 

where {/i^ : k — 0,1,2, ... , N u } is a partition of the interval 
(— 1, 1), {(j)£ : I = 0, 1, . . . , N,p} is a partition of the interval 
(0, 2ir), and the quantity P(y m , [if., tfii) A [if. A<fi is obtained 
by including only those terms on the right-hand side for which 
the direction {[i, (f>) n at time n lies in the interval [i^ < fi < 
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It is also useful to define the 



angular distribution P(y, fi) such that 
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k=0 



An At N ~ 1 

P{v m ,[i k ) A[i k = V \F{s m ,t n )\ 2 , (24) 
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where P(v m , [if.) A/i^ includes only those terms on the right- 
hand side for which the direction ([i, <p) n at time n satisfies 
fi k < fi < [i k+1 . That is, 



P{v m , [ik) A[i k 



AnAt 
Clu T 
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\F{s m ,t n )\ 2 . (25) 
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The function P(v, (i) describes the distribution of power as a 
function of the frequency v and angle 9 between the mean 
magnetic field and the average flow direction of the solar 
wind. It is the natural way to define the power spectral density 
per unit frequency and per unit angle d[x. Consistent with the 
definition of power spectral density, the total power is equal 
to the sum over all frequencies of Av m times (124) . There- 
fore, the contribution to the total power when the local mean 
field lies in a direction between (if. and [ik+i is the sum over 
m of Av m times d25l >. where the mean field depends on the 
scale s m . A different definition of the power spectral density 
is given in the next subsection. 

4.2. Power spectrum as a function of 9 

At a given scale, the number of times the local mean mag- 
netic field Bq points in a given direction r within a small 
solid-angle dfl will vary depending on the direction of the 
unit vector r. Thus, the probability distribution of directions 
of Bq on the surface of the unit sphere is not necessarily uni- 
form. By summing the energy of the fluctuations at each time 
step t n when Bq points in a given direction r, the average 
power of the fluctuations when Bq points in that direction is 
the total energy divided by the length of time Bq points in 
that direction. This is the power spectral density of the fluc- 
tuations at a given frequency when the local mean magnetic 
field points in a given direction. 

Thus, for an angle bin with direction ([i, (f) and solid angle 
d£l <C 47r, the average power per unit frequency when the 
local mean magnetic field points in that direction P(y, [i, <f) 
is defined by 



p{v m ,[i,4>) 



An 



Y, \F(s m ,t n )\ 2 , (26) 



{n: (/i,0)„ £ bin} 



where iVbin is the number of terms in the sum, that is, the 
number of times t n for which the local mean field lies within 
the given solid-angle. The definition d26b is believed to be 
equivalent to that employed by Horbury et al. (2008) and it is 
also the definition used in the present study. 

Note that the usual Fourier analysis of solar wind time se- 
ries ignores the fact that the local mean field spends differ- 
,2 ent lengths of time pointing in different directions. This un- 
' ' even distribution of directions of Bq on the surface of the 
unit sphere causes an uneven weighting of power in a stan- 
dard Fourier spectral analysis. The spectral analysis based on 
equation d26| i removes this uneven weighting by dividing by 
the number of times the local mean field points in a given di- 
rection. This is made possible by the fact that the wavelet 
transform is localized in both time and frequency and, in this 
respect, the wavelet transform is superior to the Fourier trans- 
form. 

The solar wind data listed in Table 1 was analyzed by means 
of the wavelet transform ( fT3l and the spectral decomposi- 
tion d26b . The analysis shows that for frequencies throughout 
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FIG. 4. — The number of bin counts for each angle bin at the frequency 
0.716 Hz using the data from March 2008 and January 2008 in Table 1. The 
direction 8 = is radially outward from the sun. Note that these two data 
sets are contained in opposite magnetic sectors (inward and outward). 



the inertial range and the dissipation range, frequencies from 
5 x Hz to 2 Hz in the spacecraft frame, the power spec- 
tral density in equation ( 126*1 1 is roughly independent of the az- 
imuthal angle <f) as is expected for turbulence that is cylindri- 
cally symmetric about the local mean magnetic field. There- 
fore, instead of equation ( [261 ), it is more expedient to use the 
power spectral density as a function of angle defined by 
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\F(s m ,t n )\ 2 , (27) 



where the sum on the right-hand side only includes those 
terms for which the direction (/i, <f>) n at time n satisfies Hf. < 
li < Hk+i an d N m k is the number of all such terms. The func- 
tion P{v, fi) describes the power spectral density as a function 
of the angle 8 between the local mean magnetic field and the 
(radial) mean flow direction. 

In this study, angle bins are defined by the partition 8^ — 
kir/30, k = 0, 1, 2, . . . , 30, /ifc = cos(#fc). The width of each 
angle bin is 6 degrees which is small enough to resolve the 
interesting behavior near the endpoints 8 = and 8 = n, 
but large enough to usually provide a reasonable statistical 
sample at all angles of interest. Figure |4]is a typical example 
of the distribution of bin counts seen for the data in this study. 
The orientation of the interplanetary magnetic field, outward 
or inward, is also apparent from Figure H] The maximum bin 
count coincides with the direction of the Parker spiral in each 
case. The distributions of bin counts was qualitatively and 
quantitatively similar at all frequencies from 5 x 10 -3 Hz to 
2 Hz. The number of bin counts in the parallel direction was 
typically 4%-9% of the maximum number of bin counts taken 
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FIG. 5. — The power spectral density versus frequency for angle bins cen- 
tered at 8 = 3 (bottom), 9, 15, 21,. . . , 93 degrees (t op) computed using the 
February 2008 data in Table 1 by means of equation )27t . 



over all bins. 



5. RESULTS 



The power spectra P(v, fi) computed by means of equa- 
tion d27"l ) using the February 2008 data in Table 1 are shown 
in Figure [5] The different curves in Figure [5] correspond to 
different angle bins and have been offset vertically for easier 
viewing. The minimum and maximum wavelet scales are ap- 
proximately 0.48 s and 200 s, respectively. The number of 
different scales is M = 30. The corresponding frequency 
range extends from approximately 5 x 10~ 3 Hz to 2 Hz in 
the spacecraft frame. The transition from the low-frequency 
inertial range to the high-frequency dissipation range is indi- 
cated by the change in spectral slope around 0.4 Hz, a typical 
break-frequency near the orbit of the earth at 1 AU. 

Restricting attention to the range of frequencies called the 
inertial range, v < 0.2 Hz, linear least squares fits are per- 
formed in log-log space to find the best fit power-law expo- 
nents. For a given data record, the frequency range used to 
obtain the fits is the same for all angles. But, for each data 
record in Table 1 a different frequency range is used to fit the 
data. This was necessary because of occasional outliers in 
the spectra, usually at the lowest frequencies, which can sig- 
nificantly effect the power law fits. Such outliers are partly 
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attributable to the smaller bin count at extreme angles and 
low frequencies. The frequency ranges used to fit the data 
are listed in Table 1. An example of the power-law fits are 
shown by the red line segments in Figure [3] which have been 
drawn so that they cover the precise frequency interval used 
to determine the fit. 

The dependence of the inertial range power-law exponents 
on the angle 9 are shown in Figure [6] For outward magnetic 
sectors, the power-law exponent changes from roughly 2 to 
1.6 as the angle 9 increases from to 90 degrees. The results 
are qualitatively similar in every case. For inward magnetic 
sectors, the power-law exponent changes from roughly 2 to 
1.6 as the angle 9 decreases from 180 to 90 degrees. Once 
again, the results are qualitatively similar in every case. The 
error bars for the power law exponents in Figure [6] are 99% 
confidence intervals based on linear regression analysis of the 
data on a log-log plot, they are not standard deviations. 

Near 9 — in outward sectors and 180 degrees in inward 
sectors the power-law exponents are usually more uncertain 
as a consequence of the smaller number of statistical sam- 
ples at these extreme angles. Nevertheless, the data indicate 
that the power law exponent when the mean field is parallel 
to the mean flow is approximately 2 ± 0.1. The value of the 
power-law exponent when the mean field is perpendicular to 
the mean flow is approximately 1.6 ± 0.1. In some cases the 
perpendicular power-law exponent is closer to 3/2 than 5/3, 
but the observations show significant variations from one an- 
gle bin to the next and from one data record to another so it 
is difficult to determine this value more precisely. Statistical 
studies may help answer this important question. 

The power spectral density P(v : 9) as a function of angle 
9 is shown in the upper plot in Figure [7] which shows that 
the power is peaked near the perpendicular direction and that 
the anisotropy of the fluctuations described by the ratio of 
the perpendicular to the parallel power Pi/P|| increases as 
the frequency increases through the inertial range. The iner- 
tial range observations in Figure [7j are consistent with direct 
numerical simulations of incompressible MHD turbulence 
which show that in fc-space the energy cascade is directed 
primarily perpendicular to the mean magnetic field and that 
the power anisotropy i ncreases as the wave n umber increases 
(Sheb alin et all 119831; lOughton etail Il994t iMatthaeus et all 
1996). This behavior is also incorporated into phenomeno- 
logical theories of strong incompressible MHD turbulence 
that ta ke wavevector anisotropy into account in a fundamental 
way dGoldreich & Sridha rllT995l [19971 lBoldvrevll2005l l2006t 
Lithwick et al. 2007; Beresnyak & Lazarian 2008| iChandranl 
2008; lPodesta & Bhattacharieell2009l) . 

In the lower plot in Figure [TJ the quantity P± is the average 
power in the two angle bins just below 90 degrees, bins 14 
and 15, and Pi is the value in bin 1. The ratio Pl/Pi found 
in the present study is quantitatively similar to the results of 
iBieber et al.l (119961) who found the ratio 1 .4 (mean value) at 
0.01 Hz and to the results of Horbury et al. (2008) who found 
a ratio of between 4 and 5 at 0.061 Hz at high heliographic 
latitudes near 1.4 AU. Smith et al. (2006) decomposed the 
magnetic field vector into orthogonal components perpendic- 
ular and parallel to the mean field and found that the average 
ratio of the power spectral density of the sum of the two per- 
pendicular components to the parallel component, averaged 
from 8 x 10 -3 Hz to 0.1 Hz, takes typical values around 10. 
It is not suprising that this variance anisotropy measured by 



ISmith et al. (2006) is quantitatively similar to the ratio Pi / Pu 
studied here, although the two physical quantities are differ- 
ent. Measurements similar to those of Smith et al.l (2006) 
were also reported by lLeamon et al.l (119991) who found a mean 
ratio around 10 at the high frequency end of the inertial range. 

An interesting new result is the observation that in the tur- 
bulent dissipation range the power ratio P± / P often exhibits 
a double peak structure as shown in the lower plot in Figure 
UJ What is usually called the dissipation range begins at the 
spectral break (change in slope) in Figure [5] the point where 
the spectra in Figure [5] deviate from their linear fits, which 
occurs around 0.25 Hz for the data from February 2008. For 
the same data, the lower plot in Figure [7] shows that the ratio 
P±/Pu decreases, attains a local minimum at around 0.4 Hz, 
and then increases again. What is most interesting is that it 
exhibits a prominent peak of large magnitude near 1 Hz be- 
fore decreasing significantly beyond 1 Hz. This feature has 
never been observed before. A prominent peak near 1 Hz is 
seen in all the intervals listed in Table 1. 

According to one school of thought, the turbulent en- 
ergy cascade of MHD-scale fluctuations makes a transition 
from an Alfven wave cascade to a Kinetic Alfven Wave 
(K AW) cascade at the perpendicular wavenumber k \ pi ~ 
1 (IOuataertlll998l; lOuataert & Grazinovl [19991 iLeamon et al.1 
| 1999t ICranmer & van Ballegooiienl |2003b iHowesI 120081; 
Howes et al.ll2008l;ISchekochihin et aljl2008l 120071) . Observa- 
tional evidence for such a transition i n the s olar wind is found 
in the measurements of Bale et al] (120051) which show the 
phase speed of the waves in this regime are consistent with the 
dispersion relation for KAWs. For the February 2008 data in 
the present study, the wavenumber such that k±pi ~ 1 occurs 
at the observed spacecraft frequency v sc — V sw /X± ~ 0.5 
Hz. Thus, the onset of the KAW cascade occurs near the spec- 
tral break (~ 0.3 Hz) and also near the location of the local 
minimum in the lower plot in Figure |7j The observations in 
the lower plot in Figure|7jshow that the ratio P± / P|| increases 
as the frequency increases from 0.5 Hz, reaches a peak around 
1 Hz, and then decreases rapidly near 2 Hz. The observed 
peak with the property P±/Pn ^ 1 demonstrates the perpen- 
dicular nature of the cascade in the dissipation range which is 
consistent with the existence of a kinetic Alfven wave cascade 
in this regime. 

The rapid decrease in the power ratio Px/Pl beyond 1 
Hz could be due to the dissipation of kinetic Alfven waves 
(KAWs). The linear dispersion relation for electromagnetic 
waves in a collisionless electron-proton plasma at thermal 
equilibrium contains two branches which correspond, in the 
limit k± — > 0, to the left circularly polarized ion-cyclotron 
(Alfven-cyclotron) wave and the right circularly polarized 
electron-cyclotron (magnetosonic-whistler) wave. Suppose 
that the magnetosonic-whistler can be neglected in the first 
order of approximation. For the relevant solar wind param- 
eters f3± p = 2, Vth,p = 75 km/s, and k±/k\\ — 10, the hot 
plasma dispersion relation shows that the damping of KAWs 
becomes strong, 7T ~ —1, when k±pi ~ 4. This corre- 
sponds to an observed frequency of 2 Hz where, assuming 
the wave amplitudes are sufficiently small, the kinetic Alfven 
waves are extinguished. Strong damping of KAWs near 2 Hz 
may cause the attenuation of the peak in the lower plot in 
Figure|7j However, another contributing factor may be a pop- 
ulation of parallel propagating waves as explained later in this 
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section. 

The measured plasma parameters used to make the above 
estimates for the February 2008 data in Table 1 consist of 
the average solar wind speed V^ w = 625 km/s, proton den- 
sity rip — 3 cm~ 3 , proton thermal speed v t h, P = 75 km/s, 
magnetic field magnitude B = 4.2 nT, and the Alfven speed 
Va = 49 km/s. Measurements of the electron core temper- 
atures on the Stereo spacecraft are unavailable due to instru- 
ment problems (P. Schroeder, private communication 2008). 
The ratio T v /Te = 2 is a typical value in h igh-speed wind 
dMarschl 1 99 1 at ISchwennll2006t lMarschll2006l) . however, here 



it is assumed for simplicity that T p /T e — 1. The ratio 
k±/k\\ = 10 in the last paragraph comes from the assumption 
that MHD-scale turbulence in the solar wind obeys the critical 
balance hypothesis so that ~ Sv^/va', solar wind mea- 

surements indicate that at the smallest inertial range scales 
Svk/vA ~ 1/10. The exact dispersion r elation given by equa- 
tions (10-57) and (10-66) in lStixl (119921) were used to compute 
the damping rates and Gautschi's al gorithm was use d to com- 
pute the plasma dispersion function ( Gautschi l 1 9701) . The root 
of the dispersion relation where the damping becomes strong, 
7 T ~ -1, is given by (to /VL p ) / {k\\ Pi ) = 1.5602 - 0.2375i, 
where k\\Pi = 0.4, k±pi = 4, and the proton cyclotron fre- 
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FIG. 7. — The power spectral density P(u, 8) normalized to the value in 
the parallel direction 9q = 3 degrees and plotted as a function of the angle 8 
between the measurement direction and the direction of the local mean mag- 
netic field for three different inertial range frequencies (upper plot); February 
2008 data. The anisotropy described by the ratio of the perpendicular to the 
parallel power, Pj_/Pn, increases as the frequency increases through the in- 
ertial range and a well defined peak is seen in the dissipation range at 1 Hz 
(lower plot). The red line fit is proportional to v 044 . 



quency is Q p /(27r) = 0.064 Hz. It is known that damping co- 
efficients can be sensitive to small changes in the distribution 
functions which are usually far from thermal equilibrium in 
the solar wind. Therefore, the above calculation of damping 
coefficients of KAWs in the solar wind is subject to possibly 
large errors. 

Another noteworthy feature in the lower plot in Figure [7] is 
the decrease of the ratio P± /Pn around the spectral break and 
the local minimum near 0.4 Hz. This arises from an increase 
in the power Pn measured parallel to Bq relative to the power 
P± measured perpendicular to B and can also be seen in Fig- 
ure [5] where it appears as a slight enhancement in the power 
spectrum near the spectral break in the bottom two curves. 
The probable cause is an enhanced population of nearly par- 
allel propagating waves with wavenumber knpi ~ 1. Be- 
cause the enhancement in the power spectra is only seen at 
small angles in Figure |3J 9 -C 7r/2, it is suspected that the 
waves are parallel propagating. More extensive data analysis 
is necessary to confirm this. At the present time, the source 
of the waves is unknown. It is not known if they are part 
of the turbulence, somehow generated by the turbulence, or 



whether they are independent of the turbulence. The waves 
could be generated in-situ by a kinetic instability. Coherent, 
parallel propagating waves around 0.5 Hz we re frequently de- 
tected in the solar wind by Behannon ( 1976) who tentatively 
identified both ion-cyclotron and electron cyclotron waves in 
different data se ts; similar observations have been made by 
iJian et all (120081) . Assuming the instability only exists for a 
narrow range of kn, then as kn increases beyond this range Pn 
would be expected to decrease rapidly and this may explain 
why the ratio P±/P\\ in the lower plot in Figure [7] ramps up 
rapidly to form the peak near 1 Hz. These possibilities require 
further investigation. 

Remarkably, a similar enhancement of proton density fluc- 
tu ations in the s olar wind n e ar k \\ pj ~ 1 has been reported 
bv iNeugebaueri (fl97l [19761) . iHollwegj (fl999l) has suggested 
these enhanced density fluctuations could be caused by the 
compressibility of KAWs in the case k± kn. How- 
ever, if the enhancement in the present study is caused by 
parallel propagating waves, then it cannot be explained by 
Hollweg's mechanism since, for pa rallel propagating waves, 
kj_ <C fcii . Neugeba uer et alj dl978l) have suggested that the 
enhancement could be caused by right-hand ion cyclotron 
waves driven unstable by a proton thermal anisotropy with 
Tu p > T± p . While this could also explain the wavelet obser- 
vations in the present study, other kinetic instabilities are also 
possible. The right-hand electroma gnetic (whistler) wave in- 
stability driven by a proton beam (Montgo mery et ail 1 19761) 
was dismissed by Neugebauer et al. (1978) because they be- 
lieved in this case the unstable waves exist over a frequency 
range that is wider than what is seen in observations. The 
crucial quantity, however, is not the frequency in the plasma 
frame but the wavelength because the frequency observed in 
the spacecraft frame is determined by the rapid advection of 
the waves past the spacecraft (by the solar wind) and, there- 
fore, is determined by the wavelength. Thus, by Taylor's 
hypoth esis, v sc — V sw /X. Figure 4a in Montgome ry et alj 
(1976) shows that the wavenumber range of the unstable 
waves could explain the observed density and/or magnetic 
field enhancements near knpi ~ 1. 

The rapid attenuation of the peak around 2 Hz is another 
characteristic feature of the dissipation range data in the lower 
plot in Figure|7]that could be caused by an enhanced popula- 
tion of predominantly parallel propagating waves. Evidence 
comes from the power spectra in Figure|5]which show a "flat- 
tening" of the spectra near 2 Hz or, equivalently, knpi ~ 4. 
This flattening is more pronounced at small angles but is 
present in angle bins 1 through 8 indicating that if such waves 
are present, then they may propagate at angles as large as 
8 — 45 degrees or so, even though the wave power is great- 
est for propagation near 6 = 0. It should be noted that the 
flattening of the power spectra cannot be caused by aliasing 
which would affect the spectra at all angles in the same way. 
Moreover, the underlying measurements were obtained with 
a cadence of 32 vectors per second. Further analysis of the 
data is necessary to investigate whether parallel propagating 
waves are present around 2 Hz and, if present, to characterize 
their properties which may be similar to those identified by 
Behannon (1976). 

Moving now to the inertial range, an interesting observation 
in the lower plot in Figure [7] is the red trend line with slope 
0.44 which is the best least-squares fit to the inertial range 
data between 0.025 Hz and 0.2 Hz. As discussed in the con- 
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elusions, the slope 1/2 is predicted by the phenomenological 
theory of Boldyrev (2006) and the slope 1/3 is predicted by the 
theory of Goldreich & Sridhar (1995). The best least-squares 
fit over the range from 0.01 Hz to 0.2 Hz yields the expo- 
nent 0.36 (not shown), however, it is not as good a fit as the 
exponent 0.44 over the range 0.025 Hz to 0.2 Hz. These con- 
siderations show that the scaling law is difficult to determine 
from the data. This difficulty is caused in part by the uncer- 
tainties in the ratio Px/P\\ which increase as the frequency 
decreases. These uncertainties can be reduced by using larger 
record lengths, but the record length is limited by the dura- 
tion of the high speed stream. Further study is needed of this 
important scaling law. 

6. COMPARISON TO ULYSSES RESULTS 

It is of interest to compare the results of the present study 
with those obtained by Horbury et al. (2008). Unfortunately, 
the numerical results of Horbury et al. (2008) were not made 
availible to us. Therefore, the analysis of Horbury et al. 
(2008) was repeated using the similar analysis technique de- 
veloped in the present work. 

The data consists of measurements of the solar wind mag- 
netic field by the vector helium magnetometer on the Ulysses 
spacecraft (Balogh et al. 1992). The data spans a 30 day in- 
terval in 1995, day 100 through day 129, the same interval 
analyzed by Horbury et al. (2008). The data has a nominal 
1 second cadence, however, some 1/2 s data and data gaps of 
various sizes are also present. The 1/2 s data is downsampled 
to 1 s and all data gaps are filled using linear interpolation to 
give a continuous time series with a 1 second cadence. This is 
necessary to compute convolutions using FFT techniques as 
described in section 4. The duration of the data gaps is neg- 
ligible compared to the record length and, therefore, the lin- 
early interpolated data makes a negligible contribution to the 
total power. The mean value of B = (Br, Bt, Bm) taken 
over the entire record is (1.65, —0.74, —0.03) nT. The mean 
value of B = \B\ is 3.07 nT and the r.m.s. value of B is 
5B = 2.53 nT. Thus, for the Ulysses data SB/B ~ 1. For 
each of the time intervals listed in Table 1, SB/B ~ 0.8. 

The wavelet analysis was performed in the same manner de- 
cribed in sections 4 and 5. Horbury et al. (2008) average the 
energy in each solid angle bin and then perform an azimuthal 
average with respect to (f>. In the present study, intermediate 
averages over each solid angle bin were omitted and, in accor- 
dance with equation d27l l. only averages with respect to (j> are 
performed. This has the advantage of yielding more robust 
statistics for angles in the parallel and perpendicular direc- 
tions where the bin counts are relatively small. The results of 
the Ulysses analysis are shown in Figure [8] 

Overall, the Ulysses results for the power-law indices in 
Figure [8] are qualitatively and quantitively similar to the re- 
sults for high-speed streams in the ecliptic plane shown in 
Figure [6] The power-law exponent closest to 9 = in Figure 
[8]is not as close to 2 as in the results reported by Horbury et al. 
(2008), possibly because the bin size used here is larger. The 
power-law exponents approach values near 1.6 in the perpen- 
dicular direction, notably different than the value 5/3 reported 
by Horbury et al. (2008). When intermediate averages over 
each solid angle bin are omitted and averages are performed 
only with respect to <fi, as in equation d27t . the results take the 
form of the red curve in the upper plot in Figure [8] Note that 
the red curve gives a more monotonic functional form with 
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FIG. 8. — The power law exponents as a function of the angle 6 between 
the local mean magnetic field and the mean flow direction (radial direction) 
for the 1995 Ulysses data, DOY 100-130 (upper plot). The black curve is 
obtained by averaging the power in each direction (/i, <j>) separately and then 
performing a (/(-average (upper plot). The red curve is obtained by averaging 
with respect to <j> only which yields better statistics (upper plot). The ratio 
of the perpendicular to the parallel power, P±/Pn, increases approximately 

like i^ 1 / 3 as predicted by the Goldreich & Sridhar theory (lower plot). The 
linear least-squares fit yields the slope 0.35 ± 0.04. 



a perpendicular power law exponent closer to 1.58 than 1.67. 
The precise value of the parallel and perpendicular power-law 
exponents are important for comparisons with turbulence the- 
ories. 

The ratio of power in the perpendicular direction to the 
power in the parallel direction, Px/P\\, is shown in the lower 
plot in Figure [8] Here, Pii is the power in angle bin 1, P± 
is the average power in bins 14 and 15, and the power is 
defined by equation d27| i. As described in more detail in 
the conclusions, the Goldreich & Sridhar (1995) theory pre- 
dicts Px/Pn oc v 1 ^ and the Boldyrev (2006) theory predicts 
Px/P\\ oc v 1 ! 2 . The data in the lower plot in Figure [8] yield 
the power law fit j, - 35 ± - 04 . This result is consistent with the 
Goldreich & Sridhar (1995) theory but inconsistent with the 
Boldyrev (2006) theory. These results are new and were not 
analyzed this way by Horbury et al. (2008). 

7. COMPARISON TO POWER SPECTRUM USING ALL DATA 

It is of interest to compare the results for the power law 
exponents in section 5 to the power law exponents obtained 



13 



10" 



D- 10 

s 
> 



"T — I I I 1 1 III 1 — I I I Mill 1 — I I I I Mil 1 — I I I I III 



oc=1.60 



A AAA n ooo°oo°-°°° 00000 ° 00 ° A d 



oo° 



02/2008 



6 



a=1.67 



6 -I 



uL 



uL 



_l 



_l 



10 

10" 4 10"' 10"^ 10" 1 10 U 
Spacecraft Frame Frequency (Hz) 

FIG. 9. — The compensated power spectrum v a P(v), trace spectrum, ob- 
tained using data for February 2008 in Table 1. The best fit power-law ex- 
ponent over the frequency interval from 10~ 3 Hz to 10 _1 Hz is a = 1.60 
(triangles). The exponent a = 5/3 is not a good fit (circles). The horizontal 
line is drawn for comparison. 



using all data. The power spectral density for each record 
(trace power) is computed using equation (T% which yields 
the same result as the power spectrum comp uted using stan- 
dard FFT techniques (Per cival & W alden 1993). The spectral 
exponent is obtained by a linear least squares fit over the in- 
terval from 10~ 3 Hz to 10 _1 Hz. An example of the compen- 
sated spectrum for the February 2008 data in Table 1 is shown 
in Figure|9] The power law fits for the other inervals are qual- 
itatively similar to that shown in Figure [9] A summary of the 
power law exponents for each interval are listed in Table 1 . 

The power law exponent for each interval (all data) is in ap- 
proximate agreement with the power-law exponents observed 
over the range of angles between 6 = tt/2 and the direction 
of the Parker spiral. This is to be expected because observa- 
tions in Figure [6] shown that the power law exponent is ap- 
proximately constant between 9 = tt/2 and the direction of 
the Parker spiral and the amplitude of the power spectrum is 
a maximum roughly when 9 = ir/2. Furthermore, the fact 
that the direction of the local mean magnetic field is most of- 
ten found near the direction of the Parker spiral means that 
the Fourier spectrum for all data is more heavily weighted by 
those fluctuations occuring when Bo lies near the Parker spi- 
ral direction. Unfortunately, the large uncertainties for the ex- 
ponents in Figure [6]preclude more precise comparisons of the 
power law exponents as a function of angle 9 and the power 
law exponents obtained using all data. 

8. DISCUSSION AND CONCLUSIONS 
8.1. Summary 

Wavelet analysis of solar wind data obtained at high helio- 
graphic latitudes by the Ulysses spacecraft led Horbury et al. 
(2008) to conclude that in the inertial range the power-law ex- 
ponent of solar wind magnetic field fluctuations changes from 
approximately 5/3 to approximately 2 as the angle 9 between 
the local mean magnetic field and the mean flow direction de- 
creases from 7r/2 to 0. The similarity between high-speed so- 
lar wind at high latitudes and high-speed streams in the eclip- 
tic plane around solar minimum suggests that a similar effect 



should exist for power-law exponents in high-speed streams. 
The results of this study confirm this expectation. 

Each of the nine high-speed streams analyzed in this study 
showed a decrease in the power-law exponent from 2 ± 0.1 to 
1.6 ± 0.1 as the angle 9 increases from to tt/2 for streams 
embedded in outward-directed magnetic sectors (large scale 
Bo pointing away from the sun) and as the angle 9 decreases 
from 7r to 7r/2 for streams embedded in inward-directed mag- 
netic sectors (large scale Bq pointing toward the sun). 

In the inertial range, observations show that the power in the 
perpendicular direction 9 = ir/2 is greater than the power in 
the parallel direction 9 = or 9 = tt by factors ranging from 
approximately 2 at 0.01 Hz, near the middle of the inertial 
range to a factor of roughly 10 at 0.2 Hz, the upper end of the 
MHD-scale inertial range. Thus, it may seem that the power 
spectrum of the entire data record should have approximately 
the same power-law exponent as the spectrum in the perpen- 
dicular direction 9 = tt/2. If true, then this provides a way 
to obtain more accurate estimates of the power-law exponent 
when 9 = tt/2. However, as discussed in section 4.2, be- 
cause the usual Fourier frequency spectrum is weighted by an 
uneven distribution of directions of the local mean magnetic 
field which is peaked in the direction of the Parker spiral, the 
power spectrum of the entire data record will likely have a 
power-law exponent that lies between the values observed in 
the directions 9 = tt/2 and the direction of the Parker spiral, 
and probably closer to that in the direction of the Parker spiral 
because of its heavier statistical weight. 

Interestingly, this effect could cause a spectrum with a per- 
pendicular power-law exponent near 3/2, as measured using 
the wavelet technique at 6 = tt/2, to produce a power-law ex- 
ponent that is close to 1.6 or even 5/3 as measured by Fourier 
spectral analysis of the entire record. In fact, the results in 
Table 1 and in Figure [9] show that the power law exponent 
for the total magnetic energy is closer to 1.60 than 1.67 in all 
the inervals studied here. Therefore, it is possible that Fourier 
analysis techniques and Fourier frequency spectra which have 
been a standard tool in solar wind research for five decades 
may be giving misleading results for the power law index of 
turbulent MHD-scale fluctuations. 

8.2. Comparison to previous work 

It is o f special interest to c omment on the large statistical 
study by iTessein et al.l (120091) who examined both high- and 
low-speed wind in the ecliptic plane near 1 AU and found no 
statistically significant dependence of the spectral index on 
the angle 9. The analysis technique employed by Tessein et 
al. (2009) was based on 1 hour intervals of solar wind data 
with a uniform cadence of 64 seconds. The spectral exponent 
for each 1 hour interval was determined by means of second 
order structure functions and the mean magnetic field for each 
interval, a 1 hour average, was used to bin the data according 
to the angle 9. 

Tessein et al. (2009) may have been unable to detect any 
angle dependence of the spectral index, in part, because they 
neglected the scale dependence of the local mean magnetic 
field. This is important because the dynamically relevant 
mean magnetic field for fluctuations at the 10 second scale 
is not the mean magnetic field at the 1 hour scale. Note, how- 
ever, that even if they did use a scale dependent mean mag- 
netic field, the 1 hour intervals they use probably contain in- 
sufficient data to measure the spectral index near the parallel 
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direction 6 = 0, at least for the high-speed streams studied 
here. To demonstrate this, the wavelet analysis technique was 
applied to data intervals with a duration of 2 hours, subinter- 
vals of the data listed in Table 1 . The 2 hour intervals are com- 
parable in length to the 1 hour intervals used by Tessein et al. 
(2009). It was found that the two or three angle bins closest to 
8 = usually had insufficient data to form a power spectrum. 
Therefore, when restricted to such short record lengths, the 
wavelet technique is unable to detect changes in the spectral 
exponent near 8 = 0. Even with 8 Hz data, record lengths 
of a few hours do not contain sufficient data near the paral- 
lel direction to reveal the change in spectral index. Like the 
wavelet analysis, the analysis based on second order structure 
functions should be able to detect the angle dependence of 
the spectral index provided a large enough statistical sample 
is used and provided the scale dependence of the local mean 
magnetic field is taken into account. 

There appears to be a di screpancy between t he work re- 
ported here and the work of iDasso et alJ (120051) . iDasso etaf] 
(2005) used data from the Advanced Composition Explorer 
(ACE) in the ecliptic plane near 1 AU to study correlation 
functions of solar wind fluctuations as functions of the an- 
gle 8 and concluded that in high-speed streams the power in 
the fluctuations is dominated by wavevectors nearly parallel 
to the local mean magnetic field. This result of Dasso et al. 
(2005) appears to be inconsistent with the results presented 
here which show that for each of the high-speed streams in- 
vestigated P±/Pii > 1 throughout the inertial range, that is, 
the power when the spacecraft traverses a path perpendicu- 
lar to Bq is greater than the power when the spacecraft tra- 
verses a path parallel to Bq. This implies that power in the 
3D wavevector spectrum is dominated by fluctuations with 
fcj_ 3> fen. This is also a well known featur e of MHD tur- 
bulence with a strong mean magnetic field dShebalin et alJ 
119831; lOughton et alj|1994t iMatthaeus et al.ll 19961) . It should 
be noted that the study of IDasso et al.1 (12005b encompasses 
larger inertial range scales than those considered here. How- 
ever, it is not known if this can explain the discrepancy. 

Hamilton et al.1 d2008l) studied the variance anisotropy of 
magnetic field fluctuations (5B±) 2 /(5B\\) 2 , the ratio of the 
power in the two components of the magnetic field vector per- 
pendicular to Bo to the component parallel to Bq. Their re- 
sults show that for (3 p ~ 1 the ratio (5B±/5B\\) 2 takes typical 
values between 1 and 5, both in the inertial range and dissipa- 
tion range. Moreover, in the dissipation range this ratio is typ- 
ically less than or equal to that in the inertial range indicating 
an increase in the parallel component relative to the perpen- 
dicular components. This behavior is qualitatively consistent 
with theoretical predictions for KAWs (Hollweg 1999). The 
power ratio P_\_/P\\ in the present study is different from the 
variance anisotropy investigated by Hamilton et al.l (120081) be- 
cause P± and Pu both represent the trace power (three vector 
components). The wavelet technique can easily be adapted to 
measure the powers of the magnetic field components and the 
associated variance anisotropy. This is an interesting avenue 
of investigation for future work. 

8.3. Theoretical interpretation 

The observed variations of the spectral exponent as a func- 
tion of 8 have a theoretical interpretation that sheds light on 
the 3D wavevector anisotropy of solar wind turbulence in high 
speed wind. Horbury et al. (2008) have suggested the results 
may be interpreted as evidence that solar wind turbulence is 



characterized by a 3D wavevector spectrum of the form given 
in the Goldreich & Sridhar (1995) theory of strong incom- 
pressible MHD turb ulence, a theo ry based in part on the fun- 
damental work by Higdon] (119841) . This is because the func- 
tional form of the Goldreich-Sridhar spectrum implies that the 
reduced spectra, the spectra observed in the spacecraft frame, 
have power-law exponents of 5/3 and 2 in the perpendicular 
and parallel directions, respectively; see Tessein et al. (2009) 
for a derivation of this result. However, it should be noted that 
Boldyrev's (2006) theory of strong incompressible MHD tur- 
bulence gives a similar prediction with power-law exponents 
of 3/2 and 2 in the perpendicular and parallel directions, re- 
spectively; the appropriate form of the wav evector spectrum 
in this case is discussed by Podestaj (120091) . The analysis of 
solar wind data shows that the perpendicular power-law ex- 
ponents 3/2 and 5/3 cannot be distinguished from the data 
because of experimental uncertainties. Therefore, it may be 
premature to draw inferences regarding agreement with one 
or another theory. 

Simulations of incompressible MHD turbulence suggest 
that Boldyrev's scaling applies in the presence of a strong 
mean magnetic field such that 5b <C Bq, where 5b is the 
r.m.s. level of the tur bulence and Bq is the mag nitude of the 
mean magnetic field (Mason et al. 2006} 120081) . These and 
other simulations also suggest that the scaling may revert to 
the Goldreich-Sridhar scaling when the mean magnetic field 
is weak, 5b > Bq, although Boldyrev (2006) has speculated 
that his scaling may still apply in the inertial range provided 
the inequality Sb <C Bq holds for fluctuation amplitudes 5b 
within the inertial range. For both the Ulysses data analyzed 
by Horbury et al. (2008) and the high-speed streams analyzed 
here, the r.m.s. level of the turbulence satisfies 5b ~ Bq while 
at the smallest inertial range scales 5b <C Bq. Observations of 
the power ratio P±/Pii in Figure [8] show that for the Ulysses 
data the Goldreich-Sridhar scaling gives a better fit than the 
Boldyrev scaling. Therefore, observations of this one time 
interval indicate that the Goldreich-Sridhar scaling is more 
consistent with the data even though 5b <C Bq at the small- 
est inertial range scales. The scaling law of the power ratio 
P±/P\\ for high-speed streams in the ecliptic plane is difficult 
to determine because there is insufficient data to obtain reli- 
able measurements at low frequencies. This problem can be 
resolved by using longer record lengths, however, the record 
length is limited by the lifetime of the streams. 

From a theoretical point of view, it is important to real- 
ize that neither of the above theories can be applied to so- 
lar wind turbulence, even if incompressibility is a reasonable 
first approximation for the description of large MHD-scale 
solar wind fluctuations. This is because the theories of Gol- 
dreich & Sridhar (1995) and Boldyrev (2006) both assume 
that the cross-helicity of the turbulence vanishes, a point also 
made explicit by Higdon (1984), an assumption that i s usually 
not sat isfied in the solar wind ( Marsch ~& Tul Il990t iMarschl 
1991b). A generalization of these theories to turbulence with 
nonvanishing cross-helicity is required to provide the proper 
foundation for comparisons between turbulence theory and 
solar wind observations. In particular, the theoretical form of 
the 3D wavevector spectrum of the two Elsasser spectra E + 
and E~ is needed. 

Lithwick et al. (2007) have extended the theory of Gol- 
dreich & Sridhar (1995) to imbalanced turbulence, that is, 
turbulence with nonvanishing cross-helicity. In the theory 
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of Lithwick et al. (2007) the parallel correlation lengths 
of plus and minus Elsasser wavepackets are equal, that is, 
Ajj" = AjT. However, this is inconsistent with the existence 
of different cascade times for the plus and minus Elsasser 
wavepackets as found in the theory. This inconsistency is 
easily remedied and the perpendicular energy spectra of the 

— 5 /3 

two Elsasser fields still have the same k ± scaling found 
by Lithwick et al. (2007). With this modification, the en- 
ergy cascades of the plus and minus Elsasser variables are 
both in a state of critical balance and the reduced spectra 
perpendicular and parallel to the mean magnetic field are 

predicted to scale like kj 5 ^ 3 and fc^ 2 , respectively. Like- 
wise, in the generalizations of the theory of Boldyrev ( 2006) 
to imbalanced turbulence by | Perez & Boldvrevl (120091) and 
iPodesta &~ Bhattachariee (2009J), the reduced spectra perpen- 
dicular and parallel to the mean magnetic field are predicted 

to scale like kj 3 ^ 2 and k^ 2 , respectively. The solar wind ob- 
servations of Horbury et al. (2008) and those presented here 



support these theoretical predictions although they cannot yet 

— 5 /3 

distinguish the Goldreich & Sridhar scaling k ± from the 

3/2 

Boldyrev scaling k ± based on power-law exponents alone. 



A remarkable conclusion of this study is that single space- 
craft measurements can be used to glean information about 
the wavevector anisotropy of solar wind turbulence. Whether 
the power spectrum measured in the direction perpendicular 
to the local mean magnetic field has an exponent of 5/3 or 
3/2 is an important fundamental question that remains unan- 
swered. Interestingly, the 3/2 spectrum for magnetic field 
fluctuations is not yet ruled out. 

It is a pleasure to thank Amitava Bhattacharjee, Charles 
Smith, Bernard Vasquez, Sergei Markovskii, Mario Acuna, 
David Curtis, and Jean Perez for helpful discussions. I am es- 
pecially grateful to S. Peter Gary for his interest in this work 
and for several helpful comments and suggestions. 



APPENDIX 

DERIVATION OF THE INVERSION FORMULA 

Equation (Q is derived from equation © as follows. The Fourier transform of the continuous wavelet transform (Q~|i is 



F(s,t)e~ luJt dt = \s\ 1/2 f(uj)4>*(us), 



(Al) 



where f(u>) and ip(uj) are the Fourier transforms of f(t) and ip(t), respectively. Therefore, by Parseval's theorem for Fourier 
transforms, 



F(s,t)G*(s,t) 



ds dt 



ds 



djj 



-I ^/HrHI^Ml 



(A2) 



Assuming ip(uj) = for uj < 0, the right-hand side of equation (IA2t is equal to 



ds 



diu 



2tt 



ds 



2ir 



ds 



dui 
2^ 



[/Mff» + /'(-^).9*(-^)]|^M| 2 . (A3) 



If the functions f(t) and g(t) are real valued, then f(—u>) = /*(<*>) and this becomes 



= 2Re 



ds 



duj 



-I ^JWMI^MI' 



Changing the order of integration and performing the integral with respect to s yields 

duj 



2CRe 



2tt 



f{w)g*{u) = C f(t)g(t)dt 



(A4) 



(A5) 



(write the integral on the left-hand side as an integral from — oo to +oo and then apply Parseval's theorem). Now, from ( 1A31 >. this 
is also equal to 



2 Re 



°° ds 
Isl 



du> 
2^ 



2 Re 



ds 
72" 



dtF(s,t)G*(s,t). 



Hence, 



2 Re 



ds 



dtF(s,t)G*(s,t)=C f(t)g(t)dt 



(A6) 



(A7) 



Setting g(t) = 5(t — r) yields equation © and setting g = f yields equation ©. 
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WAVELET TRANSFORM OF cos(t^T) 

Inserting f(t) — cos(ojt) into equation (Q~|) yields the continuous wavelet transform 

F{s,t) = 2- 1 / 2 ^ 1 / 4 |s| 1 / 2 [^exp(iwt) + Bexp(-iwt)] 

where 

A = exp[-(ws - lu ) 2 /2] - exp{-[(ujs) 2 +u%]/2}, 
B = exp[-(ws + ujo) 2 /2} - exp{-[(ws) 2 + w§]/2}. 
The distribution of energy in the (s, t) plane is given by the function 

2AB 



\F{s,t)\ 



|s|(A 2 + £ 2 



1 



It is easy to show that A/B — — exp(wtJos) and, therefore, 



\F(s,t)\ 



2 T 1/2 ,„,^2 , ^ 



A 2 +73 2 

cos(2cji) 



cos(2wi) 



cosh(wo;os) 



(Bl) 

(B2) 
(B3) 

(B4) 
(B5) 



The last term in the square brackets is negligible when loujqs > 3, that is, it is negligible except at very small frequencies or small 
scales. With this approximation, 



\F(s,t) 



\s\(A z + B z ). 



(B6) 



To find the scale s where this function attains its maximum, assume s > and then differentiate this expression with respect to s 
and set the result equal to zero. This yields the equation 



,'2x 2 \ cosh(3iE/2) 
4>( x ) = — - 1 cosh x - x . = 0, 

>0 / smn(x/2) 



(B7) 



where x — loloqs. This equation has exactly one root in the interval < x < oo. It has already been assumed that cosh(a;) 3> 1 
and, therefore, this may be approximated by the quadratic equation 



2a; 2 



2x - 1 = 



(B8) 



with solution x ~ luq + |. Hence, the scale s where the wavelet transform of cos(W) is a maximum is approximately 



ujujo 



This proves ( fT2l . 



(B9) 
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